Quantum decay of dark solitons 
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Unless protected by the exact integrability, solitons are subject to dissipative forces, originating 
from a thermally fluctuating background. At low enough temperatures T background fluctuations 
should be considered as being quantized which enables us to calculate flnite lifetime of the solitons 
T ~ T~'^. We also find that the coherent nature of the quantum fluctuations leads to long-range 
interactions between the solitons mediated by the superradiation. Our results are of relevance to 
current experiments with ultracold atoms, while the approach may be extended to solitons in other 
media. 
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Soliton dynamics is in a heart of multiple areas of 
physics and applied mathematics [1]. The recent resur- 
gence of interest in propagation of solitons through a dy- 
namic media [21 13| was stimulated by a number of exper- 
imental observations [H [S] of dark solitons (DS) in one 
dimensional (ID) Bose systems. The studies of the soli- 
ton dissipative dynamics [21 13 E] are probably the most 
relevant for the current experiments. Indeed the record 
lifetime of the dark soliton is ^ 200ms [5 , calling for 
understanding of the ultimate quantum limits for soliton 
persistence. 

As follows from the Mermin- Wagner theorem, the 
mean- field approach is not valid in ID even at zero tem- 
perature. This calls for an essentially quantum descrip- 
tion, first achieved by Lieb and Liniger [7] for ID bosons 
with the contact two-body interactions, solvable within 
the Bethe ansatz. As was latter shown by Kulish, Man- 
akov and Faddeev [8] one of the fundamental excitations 
within the Bethe ansatz classification (so called Lieb II 
mode jlj) essentially coincides with DS. This fact pre- 
serves existence of DS within the Lieb-Liniger model even 
on the quantum level. The underlying reason for such a 
stability is the infinite set of conservation laws, charac- 
teristic for exactly integrable models [TOj . It allows exact 
canonical transformation to collective coordinates of the 
soliton and its description as a point-like particle obeying 
quantum mechanics . At T = the Heisenberg uncer- 
tainty of the soliton position and its quantum spreading 
during the evolution impedes the imaging of dark solitons 

At finite T the soliton loses its quantum coherence 
and it was realized in Ref. ^ that the major factor 
limiting DS lifetime is the lack of the exact integrabil- 
ity. In a strongly confined (with transverse frequency 
'^-l) gas it stems from the virtual transitions to higher 
states of the transverse quantization leading to the ef- 
fective three-body contact interactions with the coupling 
constant a — — 12 ln(4/3) ^^/(fiwj,) [13]. Here g is the 
strength of delta-like two-body interaction. On the clas- 
sical level such three-body interactions lead to a cubic 
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FIG. 1; a) Dispersion relation E{P) of soliton and two- 
phonon processes leading to the dissipation. The arrows rep- 
resent absorption and emission of long wavelength phonons; 
their slope is given by sound velocity c. b) Function f{y) 
defined in Eq. (l20|. 



in density term SL = 



Y6 in GP Lagrangian [2]- 



The corresponding modified Gross-Pitaevskii (GP) equa- 
tion still possesses a one parameter family of soliton so- 
lutions, however, due to the broken classical integrabil- 
ity, such solitons are unstable against scattering on ther- 
mally excited density fluctuations. This classical mech- 
anism of soliton dissipation was shown to result in a fi- 
nite lifetime r ~ (r//i)~^, for temperatures T higher 
than chemical potential fi of ID liquid [2|. In this Letter 
we report a study of the soliton dynamics in the low- 
temperature regime T <^ fi. For such low temperatures 
the phonons are "cold": the density and phase fluctua- 
tions of the background need to be treated as quantized 
objects. Since the soliton velocity V is always smaller 
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than the speed of sound c, emission of a single phonon 
is forbidden by the energy and momentum conservation, 
i.e. by Landau criterion. The leading allowed process 
is the Raman two-phonon scattering depicted in Fig. [T^, 
where one thermal phonon is absorbed and another one 
reemitted. The rate of these processes and their physical 
consequences in the dynamics of DS constitutes our main 
results which are summarized in Eqs. (l]|4) followed by 
their derivation. 

For small soliton velocities V <^ c the characteristic 
momentum transferred to phonons q T/c is much less 
than the typical DS momentum nhno justifying the semi- 
lassical treatment and we find the equation of motion 



M*V = -nV . 



(1) 



Here M* = —4hno/c is the (negative) effective mass of 
the soliton, while the right hand side is the viscous fric- 
tion force due to phonon scattering . Due to the neg- 
ative effective mass M* the peculiarity of DS dynamics is 
that the friction force accelerates rather than decelerates 
the soliton. The accelerated DS looses its energy and 
eventually thermalizes with the phonons. The quantum 
nature of the phonons manifests itself in long-range inter- 
actions between solitons if one considers a generalization 
of Eq. ^ for a gas of DS with the set of coordinates Xi 
and velocities Vi, 
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where the temperature length is defined as Lx = hc/T 



and the function f{y) is defined in Eq. (201 and shown 
in Fig. ^p. It is characterized by /(O) = 1 and f{y) w 
—30TTye~'^^y for y ^ 1. As a result, only solitons within 
the distance Xi — Xj < Lt interact effectively with each 
other. This is clearly only possible deep in the quantum 
regime T <^ fi, where is much larger than DS size 
given by the healing length ^ = H/mc. Contrary to the 
usual potential forces, the mutual friction in Eq. ^ is 
an even function of the relative coordinate, and an odd 
function of the center of mass velocity. As a result it does 
not affect the relative motion M*{Vi — 1^2) = ^^(^2 — ^1) 
while the center of mass motion is strongly affected 
M*{Vi+V2) = -K{Vi + V2)[l+f{{Xi~X2)/LT)]. Incase 
Xi ~ X2 0.2Lt the two DS accelerate almost twice as 
fast as a single one. This is a consequence of the super- 
radiation, when phonons emitted by one DS stimulate 
coherent emission by the other. The same mechanism 
slows down the acceleration of two DS center of mass for 
Xi-X2> 0.2Lt, Fig.[l]3. 

We found for the viscosity coefficient 



''^^>~ 1215 hc^ 



(3) 



Its sublinear dependence on temperature reflects new 
regime of "cold" degenerate phonons. In the opposite 



limit T ^ /i we find linear law k ^ (a^nQ/fi.c^)(T//i) 
in agreement with Ref. |2], which is a result of the soli- 
ton interacting with "hot" classical phonons. The typical 
lifetime [TS] of DS may be estimated from Eq. ([T]) as 
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Taking ®''Rb atoms confined with transverse frequency 
~ IkHz we find q — 10"^^ J m, and —a = 



UJj_ 



J m^, which for den- 



Ulog{4:/3)g^/hiuj_ ~ 3.5 • 10" 
sity no — 10* m^^ and sound velocity c ~ 1 mm/s yields 
T ~ lOOmsx (/i/T)'*. Experimentally ^/T ^ 1 so our re- 
sults are in a reasonable agreement with the observed [5] 
soliton life-time of r '-^ 200 ms. 

We first discuss the nature of the soliton-phonon in- 
teractions. To this end we use the integrable case with- 
out the cubic in density n'^ term in the Lagrangian, and 
verify the absence of phonon emission by DS. We then 
reintroduce a and derive the dissipative force Eqs. ([T]), 
([3]). The weakly interacting Bose gas may be described 
by the Lagrangian 



L= dx 
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2m 



{n-nof , (5) 



where n{x,t), and 0(a;,t) are density and phase fields. 
The gradient of the latter is related to the superfluid 
velocity u{x, t) = dx4>{x, t) /m, where m is the mass of the 
bosons. In Eq. ^ we have subtracted the contribution of 
uniform density profile np. Variation of the Lagrangian 
([5]) with respect to the fields gives the equations of motion 



dtn = 
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(6) 
(7) 



This is the hydrodynamic form of GP equation |16j . 
Eq. ^ is the continuity equation, while Eq. ([T]) is the 
Euler equation. In the long wavelength limit, they re- 
duce to the wave equation describing low energy sound 
waves (phonons) propagating with the sound velocity 
c. The latter depends on the background density no 
and can be obtained from the thermodynamic rela- 
tion mc^ /no — d^/dno — g. The corresponding La- 
grangian written in terms of small variations of the den- 
sity p{x, t) = n{x, t) — no becomes 



Lph = / dx 



no 



2m 



{dx^f 



2no 



(8) 



The nonlinear phonon interactions term p{dx4>)'^ is re- 
tained since it is going to play an important role below. 

In addition to the phonon modes, governed by Eq. (|8|, 
equations of motion (|6]), ([t]) support a one parameter 
family of soliton solutions [T^ with the localized profiles 
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of density ns{x ~ Vt; V) and current Us{x — Vt; V) given 
by 

sin2($/2) 




cosh^ (zTOcsin(<i>/2)) 
no 



ns{z; V) 




The velocity V is related to the the phase drop across the 
soliton $ = 0(-oo) - (j){+oo) as V/c = cos($/2). One 
also finds for the momentum and energy of the soliton 

P{V) =noi<P- sin $) ; E{V) = Eq sin3($/2) , (If) 

where i?o("o) ~ 4cno/3 can be regarded as the soliton 
chemical potential. Expressions ( f f I lead to the equation 



of motion for a soliton moving with a constant velocity 
V = dvE/dyP = dE/dP. In what follows we shall focus 
on dark, $ « tt, small velocity solitons V ^ c, which have 
P{V) « TTOo + M*V and E(V) ^ Eq + M*V^/2 15J. 

We now substitute the soliton solution, rig (a; — X;V), 
and Us{x-X;V), Eqs. into Eq. ^. Here X{t) 

is an instantaneous soliton coordinate related to the ve- 
locity V by the equation of motion X = V. As a result 
one finds an effective DS Lagrangian 



Ls = P{X)X - E{X) = nnoX 



M*X' 



— En 



(12) 



It provides Feynman path integral description of the soli- 
ton as a quantum particle with the coordinate X(t) mov- 
ing in constant background with density n{x, t) — hq 
and u{x,t) — 0. Local fluctuations of density hq — > 
no + p{X, t) and background velocity u{X, t) interact 
with the soliton and modify its dynamics. In the long- 
wavelength limit the effect of phonons can be studied 
by considering small changes to the uniform background. 
The change in the background density is accounted for by 
modifying the chemical potential of the soliton, the last 
term in Eq. (|12[), by expanding it as Eo{nQ -I- p) ~ i?o + 



iVo(mc7no)p + a„„(iVomc7no)pV2. Here TVq - OEo/dfi 
is the number of particles expelled from the dark soliton. 
To obtain the coupling to the velocity field, we note that 
in the laboratory frame, where the liquid moves with the 
uniform velocity u, the fields Eqs. (|9|, ( [lO| transform as 
ns{z, V) — )• ns(z, V—u), Us{z, V) — u+Us{z, V—u). This 
transformation modifies momentum P{V) — > P(V — u) 
and energy E{V) — )■ E{V — u) + uP{V — u) of the soli- 
ton due to the uniform background flow u \17\ . Together 
with the density corrections it leads to the Lagrangian 



Ls + is-ph ^ Ls{X - u, no + p) = -7r(no + p){X ~ u) 

(13) 



M*X^ ■ M*u^ cp' 

M*Xu^ 2cp 

2 2 2no 



where p = p{X, t) and u = u{X^ t). 

Instead of tackling it directly, it is convenient to per- 
form gauge transformation of the phonon density and 



velocity fields to get rid of terms linear in these vari- 
ables. This is achieved by the following substitution 
p{x,t) — >■ p{x,t) — NoS{x — X{t)) along with u{x,t) ^• 
u{x, t) — {n/m)5{x — X{t)). One should also redefine the 
soliton coordinate [IH] - (1 mNo/M*)u{X, t) ^ X 
to account for the phonon drag. Upon this change of vari- 
ables in Lph + Ls + Lg-ph the soliton-phonon interaction 
Lagrangian acquires the form 

Ls-pi, = -^p'{X,t)-^u'{X,t), (14) 

where the soliton - two-phonon interaction vertices are 
given by 



^ ^ dp dNo_ 



r„ = mNo 1 



miVo 
M* 



(15) 



Notice that to derive Eqs. ( 14 ) and (151 we crucially used 
the phonon non-linearity mpu^/2 in the Lagrangian (|8|. 
Since we have succeeded to transform the interaction La- 
grangian to the form which does not contain terms linear 
in the phonon fields, we can disregard now the phonon 
non-linearity and treat them as the Luttinger liquid |19j , 
described by the Gaussian part of Eq. Q {pu^ term does 
not contribute to the leading temperature dependence). 

We have arrived thus at the problem of the "quan- 
tum impurity" with the mass M*, interacting with the 



Luttinger liquid through the two-phonon vertices (14 1. 
This problem was considered in Refs. |20j in close anal- 
ogy to three dimensional dynamics of impurities in liquid 
''He [H]. The leading ^ T"* contribution to the viscous 
force acting on the "impurity" corresponds to the pro- 
cess shown in Fig. [l^, where the " impurity" absorbs one 
thermal phonon, while emitting another long wavelength 
phonon to satisfy momentum and energy conservation. It 
results in equation of motion M*V = F{V,Q), where the 
friction force exerted by the liquid on moving impurity is 



e'«^gn(g,(jF).(16) 



|q|<mc 

Here H(g, lo) is the imaginary part of the Fourier trans- 
form of 6'(t)([p^(x, <), p^(0, 0)]) response function of the 
phonon gas. For small velocity, F ^ c, one finds 



Ii{q,qV) 



Sm'^c^T sinh2(cg/4r) 



V . 



(17) 



The momentum sum in Eq. (16) is limited to \q\ < mc. 
Indeed, phonons with the wavelengths shorter than DS 
size (mc)~^ practically do not interact with the latter. 

Using the relations M* = —2mNo = — 4?io/c and 
mc^/no = g we find Tp — c/uq and r„ = tiq/c. In 
view of Eq. ( 16 ) V = 0, meaning that the soliton mo- 



tion is unaffected by the interactions with the phonons. 
This is a consequence of the exact integrability of the La- 
grangian in Eq. ([5|, which protects the soliton from the 
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dissipation. This is also expected from the fact that the 
soliton configuration is an example of a reflect ionless po- 
tential, playing an important role in the classical theory 
of integrable models [TD], so it does not scatter phonons. 
Remarkably, in our approach this fact manifests itself 
through a subtle destructive interference of the phonons 
excited by the density and current vertices. 

When a small cubic in density term —a{n — uq)^ /6 is 
added to the Lagrangian ([5| , it breaks the exact integra- 
bility of the problem. Below we calculate the corrections 
to the Lagrangian p4|) and demonstrate the lack of the 



with the phonons by the Lagrangian ( 14 ). Integrating out 



exact cancelation of the prefactor in Eq. (16), leading to 



the dissipation and eventual evaporation of the soliton by 
emission of phonons. To find corrections linear in a to 



the soliton Lagrangian ( 12 ) it is sufhcient to substitute 
the bare soliton configuration, Eqs. (|9| and (10), into 
the cubic in density term. 
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dx[n,{x;V)-nof = — — ^ 
45 mc 



y2 

^2 



(18) 

Expanding in V/c <^ 1, and comparing with Eq. ( 12 ), one 
finds 6Eo = -8ano/(45mc) and dM* = -8ano/(9mc^). 
To calculate corrections to the number of particles ex- 
pelled from the soliton Nq = {nQ/mc^)dEQ/dno and its 
derivative which enters Eq. ( 15 ) it is important to take 



into account the renormalization of the sound velocity 
mc^/nQ — d^/duQ = g + auQ. We obtain the modified 
vertices 



c 

no 



2 miQ 

3 mc^ 



no 
c 



2 ariQ 
9 mc^ 



(19) 



resulting in the non-zero soliton - two-phonon coupling 
in Eq. Tp - TuC^/nl = -8ano/(9mc). 

To obtain Eqs. ([T]), ([3| we notice that at small temper- 
ature T <^ mc^ the momentum sum in Eq. ( 16 ) may be 
extended to infinity. Substituting then Eqs. ( |17[ ) and ( 19 ) 
into Eq. one finds F(V,X) = - k{T)V f {X / Lt) , 

where 



/c"* cos(2fc?//7r) 
sinh^ k 
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sinh y 



X sinh2/(3 -(- 2 sinh'^ y) — ycoshy (2 -I- cosh'^ y) . (20) 

The dependence of the friction force due to "cold" 
phonons with q < T/c <^ mc [20] is to be contrasted to 
the regime of larger temperatures T ^ fi ^ mc^ (still 
T ^ Eq) where all the phonons with q < mc contribute 
to the dissipation, their number being simply propor- 
tional to T. One finds thus k ~ {an'^/m(?y'mT in agree- 
ment with the results of Ref. |2]. The numerical coeffi- 
cient here depends on details of the momentum cutoff at 
q ^ mc and is beyond our phenomenological long wave- 
length approach. To derive the equations of motion ^ 
for the gas of DS we describe interaction of each soliton 



the Gaussian phonons and taking the semiclassical limit, 
we arrive at Eq. ^ with the force given by Eq. (20). 

In conclusion, soliton scattering on background fluctu- 
ations, leads to its decay and a finite lifetime (unless the 
model is integrable). We have shown that in the quan- 
tum regime such a lifetime is significantly longer, than 
expected from classical considerations and solitons ac- 
quire long-range mutual interactions. Our approach may 
prove useful in other areas studying solitons propagation 
in a dynamic media, such as, e.g., non-linear optics. 
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